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ABSTRACT 

We present direct measures of the ionization fractions of several sulfur ions in 
the Galactic warm ionized medium (WIM). We obtained high resolution ultraviolet 
absorption line spectroscopy of post-asymptotic giant branch stars in the the globular 
clusters Messier 3 [(/,&) = (42?2, +78?7); d = 10.2 kpc, 2 = 10.0 kpc] and Messier 5 
[(/, b) = (3?9, +46?8); d = 7.5 kpc, z = +5.3 kpc] with the Hubble Space Telescope 
and Far Ultraviolet Spectroscopic Explorer to measure, or place limits on, the column 
densities of S I, S II, S III, S IV, S VI, and H I. These clusters also house millisecond 
pulsars, whose dispersion measures give an electron column density from which we 
infer the H II column in these directions. We find fractions of S +2 in the WIM for the 
M 3 and M 5 sight lines x(S+ 2 ) = N(S +2 )/N(S) = 0.33 ± 0.07 and 0.47 ± 0.09, 
respectively, with variations perhaps related to location. With negligible quantities 
of the higher ionization states, we conclude S + and S +2 account for all of the S in 
the WIM. We extend the methodology to study the ion fractions in the warm and hot 
ionized gas of the Milky Way, including the high ions Si +3 , C +3 , N +4 , and +5 . The 
vast majority of the Galactic ionized gas is warm (T ~ 10 4 K) and photoionized (the 
WIM) or very hot (T > 4 x 10 5 K) and collisionally ionized. The common tracer of 
ionized gas beyond the Milky Way, +5 , traces < 1% of the total ionized gas mass of 
the Milky Way. 
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1. Introduction 

Since its discovery in the 1970's, the origin of the diffuse Ha emission arising from the Galaxy 
has remained something of a mystery (Haffner et al. 2009, Reynolds 1993). This emission, in 
addition to the free-free radio absorption measurements from the 1960s (see Hoyle & Ellis 1963), 
implies the existence of a diffuse distribution of free electrons outside of normal H II regions. 
This warm ionized medium (WIM), often referred to as the "Reynolds Layer" in the Milky Way, 
dominates the mass of ionized gas in the Milky Way and other galaxies. The power required to keep 
this gas ionized can be met comfortably only by the Lyman-continuum photon production by early- 
type stars (Haffner et al. 2009, Reynolds 1993). However, the ionization of this material by OB 
stars is troublesome: the large scale height of the material (~ 1 kpc or more; Haffner et al. 1999, 
Gomez et al. 2001, Gaensler et al. 2008) requires that the photons responsible for ionizing the 
WIM travel hundreds of parsecs from their point of origin. The large cross-section for absorption 
of such photons by neutral hydrogen (a ~ 6.3 x 10 -18 cm 2 ) naively implies such photons could 
only travel very small distances, of order 0.1 pc assuming a typical interstellar density (~ 1 cm -3 ). 
Various groups have presented models for the ionization of the WIM which rely on photons from 
hot stars (Mathis 2000; Sembach et al. 2000; Domgorgen & Mathis 1994), photon emission from 
cooling hot gas (Slavin et al. 2000), heating from magnetic reconnection (Hoffman et al. 2012), 
and even photons from decaying neutrinos (Sciama 1998). 

The hot star ionization models, and to a lesser extent the models invoking cooling hot gas, 
must assume that the geometry of the interstellar medium (ISM) allows for the propagation of 
ionizing photons over the large distances required (Miller & Cox 1993, Dove & Shull 1994, Dove 
et al. 2000, Wood & Mathis 2004, Wood et al. 2010). This includes the concept of density- 
bounded or "leaky" H II regions (Haffner et al. 2009). In particular, various groups have noted the 
potential importance of the complicated geometries that may exist in a turbulent, supernova-driven 
ISM (Miller & Cox 1993, Wood & Mathis 2004), or in superbubbles, where a significant amount 
of gas has been evacuated about OB associations as a result of correlated supernova explosions 
(Wood et al. 2010, Dove & Shull 1994, Dove et al. 2000). Coupled with "leaky" H II regions, 
these geometric considerations may go a long way toward explaining the propagation of ionizing 
photons from their sources to the WIM gas. However, there is not yet a complete model of the 
ionization of the WIM, and the predicted ionizing spectrum depends on the assumed sources and 
the degree to which the photons may be processed through intervening H, thereby modifying the 
source spectrum (Haffner et al. 2009). 

Very sensitive observations of forbidden metal emission lines acquired over the last decade 
with the Wisconsin H-alpha Mapper (WHAM; Reynolds et al. 1998b; Tufte 1997, Haffner 1999, 
Madsen et al. 2006) have provided data on the ionization and temperature structure of the WIM 
(Reynolds et al. 1998a, Haffner et al. 1999, Reynolds et al. 2001, Hausen et al. 2002, Madsen 
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et al. 2006). Compared with classical H II regions the WIM appears to be warmer (Madsen et 
al. 2006). Metals and helium in the WIM are generally less highly ionized than in H II regions, 
as well. In particular, the weakness of forbidden [O III] and He I recombination emission implies 
relatively low ionization fractions of +2 and He + in the WIM compared with classical H II 
regions (Reynolds & Tufte 1995, Madsen et al. 2006). Haffner et al. (1999) demonstrated that 
the ratio of forbidden metal line strengths [S II]/[N II] is sensitive to the ionization fraction of 
S + : x(S + ) = n(S + )/n(S). The typical values implied for this ionization fraction in the WIM are 
x(S + ) « 0.3 - 0.7 (Haffner et al. 1999, Madsen et al. 2006). This is higher than seen in O star H II 
regions where S +2 is more abundant than S + . The emission line constraints on WIM ionization 
imply a lower ionization parameter for the WIM compared to classical H II regions and a distinct 
spectrum, e.g., due to spectral processing of O star radiation escaping from H II regions and/or 
contributions from other sources such as cooling radiation (e.g., Slavin et al. 2000) or hot, evolved 
stars such as white dwarfs (Bregman & Harrington 1986, Haffner et al. 2009, Flores-Fajardo et al. 
2011). 

While emission line diagnostics have provided some quantitative measures of the WIM ion- 
ization, their temperature dependence adds a layer of complexity to understanding the ionization 
of the medium. Furthermore, because the excitation of these lines relies on collisions with warm 
electrons, their intensities are weighted by « n 2 e , where n e is the electron density, implying they 
may be strongly weighted toward higher density regions. In this paper we present a complemen- 
tary approach to probing the ionization state of the WIM. Following Howk et al. (2006) we use 
ultraviolet (UV) absorption line measures of the multiphase ISM along sight lines to globular clus- 
ters that also contain pulsars. The UV observations of post-asymptotic giant branch (PAGB) stars 
in the background clusters from the Hubble Space Telescope (HST) and Far Ultraviolet Spectro- 
scopic Explorer (FUSE) provide measures of metal ion column densities in these directions, while 
radio observations of pulsar dispersion measures (Hessels et al. 2007) provide determinations of 
the electron column densities. Taken together we can provide estimates of the ionization fractions 
of S°, S + , S +2 , S +3 , and S +5 in the WIM (as well as the total gas phase abundances; see Howk 
et al. 2006). We apply this technique to two extended sight lines through the Galactic WIM, one 
probing the first 10 kpc above the disk in the very local region about the sun (Howk et al. 2003, 
2006) and the other probing the first 5 kpc of the disk, but some 5 kpc projected radial distance 
toward Galactic Center (Zech et al. 2008). 

The structure of our paper is as follows. In $2] we discuss the methodology underlying our 
technique, which builds upon the discussion of elemental abundance determinations in Howk et al. 
(2006). In £|3]we discuss the UV absorption line observations, their reduction, and our assessment 
of the hydrogen and metal column densities along the two sight lines in this work. In £j4]we present 
the results of our analysis of these two sight lines and compare these with previous constraints on 
WIM physical conditions, while in $5]we discuss the implications of our results and compare them 
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with theoretical models. We give a summary of our results in $6] 



2. Methodology 

To provide a direct quantitative measure of the ionization fractions of metal ions in the WIM, 
we rely on unique sight lines that provide measures of both metal ions and neutral hydrogen from 
UV absorption lines as well as ionized hydrogen from pulsar dispersion measures (DMs). A com- 
parison of the metal ions with the H II column density derived from the pulsar dispersion measures 
gives the ionization fractions if the abundance of the metal is known. 

We write the abundance of a metal X with respect to H: 

J2N(X*) 

A ^ = N(H I) + iV(H II) + 2iV(H 2 ) ' (1) 

where N(X^) is the column density of the jth ionization stage of X and the sum is nominally 
over all ionization states. Thus, the numerator represents the total column density of the metal 
X, while the denominator, with N(E) = N(U I) + N(H II) + 2iV(H 2 ), is the total hydrogen 
column density. In most interstellar absorption line studies, A(X) is estimated by comparing the 
dominant ionization state of the element X in the warm neutral medium (WNM) with the column 
density of H I. Such studies assume that the neglect of terms in the sums in both the numerator 
and denominator of Equation \T\ has a negligible effect, which is not necessarily true (Sembach 
et al. 2000), or they attempt to make ionization corrections on the basis of models (e.g., Howk 
et al. 1999, Howk et al. 2003). Howk et al. (2006) showed that all of the significant terms 
can be accounted for along sight lines to globular clusters with UV bright stars and radio pulsars. 
Ultraviolet absorption lines provide measurements of all of the important ionic states of S and 
sometimes other metals in the WIM and WNM (§0. Those UV measurements combined with the 
pulsar dispersion measures provide information on all of the states of hydrogen. 

Once the abundance is known, the comparison of the column of a metal ion X^ that arises 
in the WIM (i.e., with no contribution from the the WNM) with the WIM H II column density 
provides a measure of the ionization fraction of X J in the WIM, x(X J ) = N(X^)/N(X). Thus, 
we can compare the ion S +2 with the hydrogen reference for the WIM, iV(H U)wim' 



c(S +2 ) N(S III) f V(H \\) WIM 

x(H+) ~ N(S) W im [ N(K) WIM 



N(S III) 
iV(H \\) WIM 



■A{S)~\ (2) 
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giving a measure of the ionization fraction of S +2 to that of H + . Here N(S III) does not require 
the WIM subscript as it arises only within the WIM, with an energy for production that precludes 
it arising in the WNM. One may write the ionization fractions for S +3 or other WIM ions in a 
similar way. We can isolate x(S +2 ) with a value of x(H + ) = 0.95 ± 0.05 based on observations 
of the weakness of [O I] emission from the Galactic WIM (Reynolds 1989, Reynolds et al. 1998a, 
Haffner et al. 2009). 

The determination of the H II column in the WIM, iV(H \t) WIM , requires two corrections to 
the electron column density, N(e~), provided by the DMs (Howk et al. 2006). The first corrects for 
the contribution to N(e~) from hot ionized gas since we cannot measure the S ions found in this hot 
gas (e.g., S VII and higher). This correction is discussed with the pulsar DMs below in §2.11 The 
other is the correction for the e contributed to the DM from ionized He. We use the helium ion- 
ization correction factor, r], discussed in detail in Howk et al. (2006) such that iV(H II) = r)N(e~). 
In this work we adopt t]wim = 0.98 ± 0.01, the "minimum helium ionization" case discussed in 
Howk et al., which assumes the ionization fraction of He +2 in the WIM is minimal. The choice 
of the minimum or maximum helium ionization correction represents a (small) systematic uncer- 
tainty. However, the maximum helium ionization case ultimately produces an inconsistency: we 
show below (even using the maximum case) that the ionization fraction of He +2 must in fact be 
very small (SI©. 



2.1. Correcting for the Hot ISM 

We estimate HIM contribution to N(e~) following the general outline of Howk et al. (2006). 
We assume the HIM electrons arise in two components, one from T ps 10 5 — 10 6 K gas traced by 
O VI the other from T> 10 6 K gas traced by X-ray absorbing/emitting gas. Thus 

N(e-) HIM = r H l iM [^(H H) 5 + N(H II) 6 ] , (3) 

where iV(H E) 5 and iV(H E) 6 are the H II columns associated with the T ~ 10 5 ~ 6 K and > 10 6 K 
gas, respectively. The helium correction factor is r] HIM . In this coronal gas, we assume all of the 
helium is in the form of He +2 , giving tjhim = 0.83 (Howk et al. 2006). 

We estimate the first term in brackets from Equation [3] from the O VI column density along 
the line of sight (Wakker et al. 2003; Savage et al. 2003). We write 

where iV(0 VI) is the measured O VI column density derived directly from FUSE data (see §3.21) . 
A(0) = (4.90 ± 0.6) x 10 -4 (Asplund et al. 2009) is the assumed gas-phase oxygen abundance, 
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and a;(0 +5 ) = 0.2 (Sutherland & Dopita 1993, Savage et al. 2003) is the ionization fraction of 
+5 in this phase of the gas (see discussion in Howk et al. 2006). This is the maximum value 
found in most models, and some of the gas probed by O VI could have an ionization fraction 
lower by a factor of a few, depending on the conditions probed by the gas. This lower-temperature 
component is the smallest contributor to the whole, representing < 25% of the total HIM column 
(in our calculations and We adopt a 50% uncertainty on the implied column given the large 
uncertainties involved. 

The T > 10 6 K gas, represented by the second term in brackets from Equation [31 has no 
direct probe along the specific sight lines studied in this work. Instead, we adopt a mean Galactic 
distribution for this gas constructed to match X-ray absorption measurements of O VII, O VIII, and 
Ne IX (Yao & Wang 2005; Yao et al. 2009). The absorption line measurements are made toward 
both Galactic and extragalactic objects. The refined model presented in Yao et al. (2009) assumes 
a thick disk structure described by a single exponential disk, Uh{z) = Uh{0) exp(— \z\/h z ), with 
a mid-plane density n#(0) = 1.4 x 10~ 3 cm~ 3 and scale height h z = 2.8 kpc. We integrate the z 
distribution of hot gas densities to give 



with the upper integration limit being the ^-height of the star. This estimate of the HIM column 
density has many simplifying assumptions whose validities are difficult to assess. As a result, we 
adopt a 50% uncertainty in our HIM column density estimates. 



In this work we apply the approach outlined above to study the ionization fractions of S , 
S + , S +2 , S +3 , and S +5 in the WEVI for two directions through the Galactic WEVI, those toward the 
globular clusters M 3 [(/,&) = (42?2, +78?7); d = 10.2 kpc, z = 10.0 kpc] and M 5 [(/,&) = 
(3?9, +46?8); d = 7.5 kpc, z = +5.3 kpc]. In what follows, we describe the UV data and analysis 
used to derive the metal and neutral hydrogen column densities along these sight lines (§3.1 through 
§3.3). We also describe the pulsar dispersion measures and their uncertainties (§3.4). 



In this work we make use of UV spectra from the Space Telescope Imaging Spectrograph 
(STIS) on board HST and from FUSE to study the metal ion and H I column densities toward 
M 3 and M 5. The sight line to M 3 is probed toward the PAGB star von Zeipel 1 128 (vZ 1 128); 




(5) 



3. Observations, Reductions, and Measurements 



3.1. UV Observations 
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the observations, reductions, and measurements for this sight line have been discussed in Howk et 
al. (2003, 2006). The ISM toward M 5 is probed along the sight line to the PAGB star ZNG 1; 
the observations and their reductions have been described by Zech et al. (2008). The latter work 
specifically studied the high velocity gas along the M 5-ZNG 1 direction. We present the first 
measurements of the low velocity gas along this sight line here. 

The STIsQ observations used in this work all employed the E140M grating to cover the spec- 
tral range ~1 150 to 1710 A at a resolution R = A/AA ~ 45,800. This provides measures of 
the ISM absorption lines at an equivalent velocity resolution Av ~ 6.5 km s _1 (FWHM). The 
M 3-vZ 1 128 observations (program ID 9150; PI Howk) were obtained with the 0"2 x 0"06 aper- 
ture. The typical signal-to-noise ratios are in the range ~ 20-40 per resolution element. The 
M 5-ZNG 1 observations (program ID 9410; PI Howk) were obtained with the 0"2 x 0'.'2 aperture. 
The two datasets have slightly different line spread function (LSF) shapes, notably differing in the 
strength of the extended wings of the LSF, due to the different apertures used. The STIS M 5- 
ZNG 1 data have signal-to-noise ratios > 25 per resolution element in the regions of interest for 
this work (Zech et al. 2008). 

The FUSE observations used here cover the spectral range 905 to 1185 A at a resolution 
R ~ 15, 000 giving a velocity resolution At> ~ 20 km s _1 (FWHM). The data were all taken 
through the LWRS 30" x 30" apertures, and we have combined observations acquired as part of 
several programs (Howk et al. 2003, Zech et al. 2008). These data have signal-to-noise ratios > 20 
per resolution element in the regions of interest. While the wavelength calibration of the STIS data 
is excellent, the FUSE absolute (and to some extent relative) wavelength calibrations are not as 
well constrained. As discussed in the earlier source papers (Zech et al. 2008, Howk et al. 2006), 
we have bootstrapped the FUSE wavelength scale to match that of the STIS data using both stellar 
and interstellar lines to determine the alignment. 



3.2. Metal Ion Column Densities 

We use the STIS and FUSE observations to derive column densities and limits for several 
metal species toward M 5-ZNG 1, with the relevant results given in Tabled] We focus on the ions 
of S, since we have access to absorption from ions that probe the WIM directly (S III and S IV). 
In addition, we give measurements of the "high ions" Si IV, C IV, N V, and O VI along both sight 
lines. The O VI is required for estimating the HIM H II column density, but we will eventually 
use these high ions to understand the ionization fractions of both the WIM and HIM in the Galaxy 



'Proffitt et al. (2002) discuss the STIS instrument characteristics in detail. 
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To determine the metal ion columns we fit the stellar continuum in the regions surround- 
ing metal absorption lines using low order Legendre polynomials. We directly integrate the line 
absorption profiles and apparent optical depths to determine the equivalent widths and apparent 
column densities, iV a , following Sembach & Savage (1992). In all cases we adopt central wave- 
lengths and oscillator strengths from Morton (2003). To calculate the limiting equivalent width 
and apparent column density of S I, we assume an intrinsic width of ~ 30 km s _1 (FWHM), the 
value derived from a single Gaussian profile fit to the S II 1250.584 A transition. All limits given 
in this work are 3a. 

Deriving the column densities of S II and S III toward M 5-ZNG 1 required special care. The 
absorption profiles for all of the transitions from these ions are shown in Figured The integrated 
apparent column densities derived from the weaker S II transitions at 1250.584 and 1253.811 A, 
which have /-values that are different by a factor of two, differ by ~ 0.07 dex, with the weaker 
transition giving a higher N a . The strongest of the S II lines observed by STIS, at 1259.519 A and 
having an /-value three times higher than the weakest, gives a significantly lower value yet for 
N a . This progression of decreasing iV a with increasing A/ suggests these lines contain unresolved 
saturated structure (Savage & Sembach 1991). 

Figure 2 compares the apparent column density profiles as a function of velocity, N a (v), for 
S II A1250, 1253, and S III A1190 (all from STIS). The weakest S II transition (at 1250.584 A) is 
shown in each of the two panels as the thin black histogram. The discrepancies between the two 
S II transitions in the top panel are consistent with expectations in the classical case of unresolved 
saturation. Since there appears to be only moderate saturation in the S II, following Savage & 
Sembach (1991) we correct the column density for S II A 1250 based on a comparison with the 1253 
transition and adopt an uncertainty appropriate for the uncertainties in the method (the integrated 
uncertainties in N a (v) being minimal). This gives our final adopted value of log iV(S II) = 15.59 ± 
0.10. This amounts to an +0.07 dex correction for saturation compared with the value derived from 
SUA 1250. 

The S III transitions toward M 5-ZNG 1 are both contaminated to varying degrees. S III 1012 
is contaminated by adjacent H 2 transitions, while the S III transition at 1 190 A is adjacent to and 
slightly contaminated by Si II at 1190.416 A (a velocity offset of +52.4 km s _1 with respect to 
S III). We deal with these contaminations in two separate ways. We correct S III Al 190 for the 
presence of underlying absorption from Si II A 1190 using an apparent optical depth profile for 
the Si II derived from the Si II transitions at 1193.290 and 1526.707 A. In each case we scale 
the relative A/-values of the transitions involved in order to estimate the contribution from the 
1 190.416 A transition of Si II as a function of velocity. Using these two comparison lines, which 
differ in strength by a factor of 3.4, leads to the same value for the integrated N a (v) profile of S III 
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A1190: logiV a (S III) = 14.66 ± 0.04. 

We assess the strength of the H 2 lines contaminating S III A 101 2 following Howket al. (2002), 
who discuss this issue for contamination of O VI absorption. We assess the strength of the H 2 lines 
contaminating S III A1012 using lines of similar /-values arising from the same rotational states. 
In this case, the derived column density is sensitive to the assumptions about this correction and 
a somewhat uncertain continuum placement. We derive column densities from S III A1012 that 
are significantly lower than those derived from the STIS-observed S III A1190: logiV a (S III) ~ 
14.53 ± 0.05 to 14.59 ± 0.04 (statistical errors only), depending on the assumptions about the 
continuum fit, details of the H 2 contamination correction, and resolution of the FUSE channel 
from which the data are adopted. The smaller columns derived from this weaker line of S III are 
smaller than those for the stronger line at 1 190 A, suggestive of some unresolved saturation in at 
least the FUSE data. However, given the uncertainties in the measurements of S III A 101 2, it is 
difficult to directly compare the two values as we did for S II above. 

On the other hand, the peak apparent optical depth of S III 1 190 as observed with the higher- 
resolution STIS instrument is only ~ 60% that of S II 1250 (including the underlying absorption 
from Si II 1 190). Even though they do not trace exactly the same gas, it is unlikely the S III 1 190 
transition is saturated to the extent of the S II 1250 transition. This, with the difference in derived 
N a (v) values between the two transitions of S III, limits the required saturation correction for the 
S III 1 190 apparent column density to be < +0.07 dex, the value adopted for S II. The only other 
comparable probes we have of the low-velocity gas are the Si IV 1393, 1402 A lines (Zech et al. 
2008), which show little to no apparent unresolved saturation for optical depths > 2x that of S III. 
However, while their shape is similar to that seen in S III, they are offset by ~ —4 km s _1 and 
trace a somewhat different mixture of gas than S 1110 Profile fitting does not ultimately provide a 
firm assessment of the possible saturation, as there is little information on the intrinsic shape of 
the S III profile from other transitions, and the H 2 contamination and unknown LSF that plague 
the FUSE observations of the 1012 A transition make that approach non-unique in our tests. We 
proceed by adopting a saturation correction of +0.035 dex and adding a systematic uncertainty 
of ±0.035 linearly to the statistical uncertainties. We adopt this uncertainty under the assumption 
that the correction should be < 0.07 dex at the 2a level. Our adopted column density is thus 
logiV(S III) = 14.70 + 0.06. 

Several other metal species are listed in Table [fl which represent more straightforward mea- 
surements. FUSE observations place stringent limits on the S IV column densities through the 
1062.664 A transition for both sight lines, while we are also able to place limits on S VI toward 



2 This appears not to be due to a wavelength calibration uncertainty, as the Si IV high velocity cloud absorption 
seen in both transitions is consistent at velocities consistent with those seen in Si II 1 190. 
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M 3-vZ 1 128 using the doublet at 933.387 and 944.523 A. For the M 5-ZNG 1 sight line, we are 
not able to place meaningful limits on S VI given the strong contamination from H 2 in the FUV. 
Given the column of S IV is well below that of S III toward M 5-ZNG 1, we assume the S VI 
column would be minimal toward this star as it is toward M 3-vZ 1 128. Limits to S I A1295.653 
absorption from STIS show its column density to be negligible along both sight lines. 

3.3. Neutral and Molecular Hydrogen Column Densities 

We derive the interstellar H I column density of these two sight lines by fitting the damping 
wings of the Lyman-a profile observed by STIS. For the sight line to M 3-vZ 1128 we adopt the 
Lyman-a derived H I column density from Howk et al. (2006). Here we determine the H I column 
density toward M 5-ZNG 1 in an identical manner. 

The neutral hydrogen column in the direction of M 5 can be estimated from H I 21 -cm 
emission observations. As noted in Zech et al. (2008), the column density derived from the 
publicly-available Leiden- Argentine-Bonn (LAB) Survey (Kalberla et al. 2005) is log iV(H I) = 
20.56 ± 0.17. This value is averaged over a 36' beam approximately centered on M 5. However, 
small-scale structure within this large beam can cause this average value to be different than that 
appropriate for comparison with the pencil-beam measurements of metal ions toward M 5-ZNG 1 
itself. We adopt an error following the recommendations of Wakker et al. (2001) for the use of 
large beam H I 21 -cm observations for absorption line studies. 

Figure 3 shows the STIS spectrum of the Lyman-a absorption line toward M 5-ZNG 1. The 
broad Lyman-a absorption is a combination of stellar atmospheric and foreground ISM absorp- 
tion. We remove the contaminating stellar absorption by using the stellar atmosphere model de- 
scribed in Zech et al. (2008) to normalize the data during our fitting procedure. The model atmo- 
sphere was calculated with TLUSTY (Hubeny & Lanz 1995) assuming an effective temperature 
T e ff = 45, 000 K, with log g = 4.48, and an atmosphere consisting of 99% He by number (W.V. 
Dixon, 2007, private communication) with a mix of abundances specified in Zech et al. (2008). 
The adopted stellar continuum is shown as the blue line in Figure 3. We use a second-order 
Legendre polynomial to match the model atmosphere to the STIS spectrum. Ideally the model 
atmosphere should provide a very good estimate of the distribution of flux with wavelength. In 
practice the models deviate slightly from the spectral energy distribution in the data. Small dis- 
crepancies between the model and observations are due to large-scale calibration and small-scale 
order combination uncertainties in the data and uncertainties in the large-scale flux distribution 
of the model itself. The polynomial parameters are treated as free parameters during the fitting 
process and contribute appropriately to the error budget. 
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The best-fit interstellar H I column density from our modeling of Lyman-a is log iV(H I) = 
20.47 ± 0.02. The best-fit profile is shown as the red line in Figure 3. This error estimate is 
dominated by uncertainties the adopted stellar properties and their affect on the stellar models, 
which we explore in the fitting process (c.f., Sonneborn et al. 2002). 

The molecular hydrogen content for both sight lines is completely negligible compared with 
the H I and H II columns. The sight line to M 3-vZ 1128 shows no detectable H 2 absorption 
(Howk et al. 2003), with log iV(H 2 ) < 14.35 (3cr) summed over J < 3. This implies a molecular 
hydrogen fraction log/ (H 2 ) = log 2iV(H 2 ) /[N(H I) + 2iV(H 2 )] < -5.3. Toward M 5-ZNG 1, 
our FUSE data show a forest of moderate strength H 2 transitions. However, none show damping 
wings, which severely limits the column densities. Deriving the precise column density is not 
crucial for this sight line, and we use the lack of damping wings and other considerations to limit 
the column: logiV(H 2 ) < 18.0 (3a) integrated over J < 5 states. This yields log/(H 2 ) < —2.4. 
The H 2 contribution to the total hydrogen column along both sight lines is negligible. 



3.4. Radio Pulsar Dispersion Measures and the Ionized Hydrogen Column Densities 

The electron column densities toward M 3 and M 5 are derived from the dispersion measures 
toward millisecond pulsars in each cluster. Following Howk et al. (2006) we average the dispersion 
measure for three M 3 pulsars (M 3A, M 3B, and M 3D), giving (DM) M 3 = 26.33 ± 0.15 pc cm~ 3 
(standard deviation). We do not include the unconfirmed pulsar M 3C in this average, although 
it gives a consistent DM (Hessels et al. 2007). The uncertainties in these DM measurements are 
very small, typically ~ 0.1 pc cm -3 or better. This average dispersion measure is equivalent to an 
electron column density of log N(e~) = 19.91 ± 0.01 in typical units (cm -2 ). For M 5 we average 
the DM for the five known millisecond pulsars (Hessels et al. 2007, Freire et al. 2008) giving 
(DM) M 5 = 29.5 ± 0.3 pc cm~ 3 (standard deviation) or log N(e~) = 19.96 ± 0.01. These are the 
total electron columns, including contributions from the WIM and the HIM along these sight lines. 

The application of Equation [3] to the sight lines in this work yields log N(e~)niM ~ 19.08 
and 19.00 cm" 2 for the M 3-vZ 1 128 and M 5-ZNG 1 sight lines, respectively (see Table©. We 
assume O VI column densities from Table [T] The result for M 3-vZ 1128 is slightly different 
than that reported in Howk et al. (2006) given the different approach to assessing the contribution 
from the hottest HIM component, iV(H II) 6 , but the difference is small enough that it produces 
a negligible change in the final results. For example, the A(S) reported here is different by only 
~ 0.01 dex compared with the earlier value. Adopting the Howk et al. methodology for the 
M 5-ZNG 1 sight line, however, would produce results that differ by ~ 0.1 dex given the greater 
importance of lower ^-height gas along that sight line. 
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4. Results 

The final results of our WIM analysis are presented in Tables [2] and [H Table[2]gives the phys- 
ical parameters derived for each of the sight lines, including the hot gas columns, WIM fractions, 
and A(S). Table [3] presents the final ionization fractions for S°, S + (see below), S +2 , S +3 , and S +5 
in the WIM for the two sight lines. To derive the ionization fractions, we adopt the S abundance 
toward M 3, which is much better determined than that toward M 5 and equivalent to the solar sys- 
tem meteoritic abundance. Sulfur does not seem to be depleted by large amounts into dust grains 
(again evidenced by the solar-like abundance toward M 3), so it is reasonable to assume A(S) is 
roughly constant in the solar neighborhood. While the M 5 sight line is toward the inner Galaxy, 
the majority of the absorption likely occurs within the first z < 1 kpc given the scale height of the 
WIM, and thus relatively close to the Sun. 

The two sight lines in this study may trace different conditions in the WIM. The high-latitude 
cluster M 3 is close to the Galactic north pole and probes the WIM in a column above the sun. On 
the other hand, M 5 lies toward the inner Galaxy at b ~ +45°, probing WIM gas somewhat interior 
to the solar circle. As shown in Table [2l the mean fractions of ionized gas along the sight lines 
differ by a factor of ~ 2. The sight line to M 3 has a total ionized gas fraction iV(H II) /N(R) = 
0.45 ± 0.090 whereas that toward M 5 is only 0.19 ± 0.03. This is likely due in large part to the 
differing paths through low-,2 gas associated with the denser, more neutral thin disk of the Galaxy. 
The M 5 sight line probes a proportionally-larger contribution from the low-2; gas. 

Within the warm ionized gas, we find ionization fractions x(S +2 ) = 0.33 ± 0.07 for M 3 
versus 0.47 ± 0.09 for M 5 (Table[3])Q Here we have transformed from x(S +2 )/x(H + ) assuming 
a hydrogen ionization fraction x(H + ) = 0.95 ± 0.05 to represent the WIM, which generally has 
x(H+) > 0.9 (Reynolds 1989, Reynolds et al. 1998a, Hausen et al. 2002, Haffner et al. 2009). 
These ionization fractions are column density- weighted averages over the entire sight lines. We 
do not detect S IV absorption along either sight line, and this absence limits the ionization fraction 
x(S +3 ) < 0.04 (3a) for both. The ionization fractions of S° in the WIM are very low, as expected: 
a?(S°) < 0.004 (3(7) along both sight lines. 

The ionization fraction of S + in the WIM cannot be directly measured through the column 
density of S II, since that ion contains significant contributions both from the WNM and WIM 
(with the WNM dominating the column). However, none of the ionization states higher than S +2 
contribute significantly to the total, with very stringent 3a limits to the ionization fractions of S +3 



The values quoted here for the M 3 sight line differ slightly from those of Howk et al. (2006) due to the difference 
in the treatment of the HIM contribution. 

4 Utilizing the directly measured value for ^4(S) toward M 5 yields x(S +2 ) = 0.67 ± 0.19. 
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and S +5 Jj Given no higher ionization states are present, the ionization fractions of S + and S +2 must 
sum to nearly unity, i.e., a;(S + ) + a;(S +2 ) ~ 1. We use this to estimate the ionization fractions of 
S+ within the WIM along these sight lines: x(S + ) = 0.67 ± 0.07 and 0.53 ± 0.09 for M 3 and M 5, 
respectively (allowing for a la contribution from S +3 in the error budget). 

The distribution of S ion fractions in the WIM as a function of ionization energy is shown in 
the top panel of Figure|4]for S° through S +5 (excluding S +4 ), with values for the M 3 and M 5 sight 
lines shown in black and red, respectively. The horizontal bars for each ion span the creation to 
destruction energies for each ionization state. On the whole the results from the two sight lines are 
in good agreement. There is a strong peak in the ion fractions for energies between ~ 10 and 35 
eV, with very little at higher energies. On the whole the two sight lines give a consistent picture of 
the ionization distribution. The M 5 sight line may slightly favor S +2 compared with that to M 3, 
but any variations are at less than 2a. 

Independent information on the metal ion fractions in the WIM can be derived from emission 
line observations. Results derived from such studies are compared with ours in Figure |5] Haffner 
et al. (2009) give a recent review of WIM emission line observations and physics (both for the 
Galactic WIM and the diffuse ionized gas in external galaxies). They combine the observations 
of Haffner et al. (1999) and Madsen et al. (2006), who used WHAM to observe forbidden metal 
line emission from the diffuse Galactic WIM. Their observations of the ratio [S II]/Ho; provide a 
measure of x(S + ) when the temperature can be estimated with observations of [N II]/Ho; (Haffner 
et al. 1999). As summarized in Haffner et al. (2009), the WHAM results give x(S + ) ~ 0.3 to 0.7, 
with the majority of the sight lines seemingly above 0.5. The results for the two sight lines studied 
here, giving x(S + ) ~ 0.5 and 0.7, are completely consistent with the WHAM estimates (Figure[5]). 

The lack of S IV and higher ion absorption is also consistent with WHAM observations of 
forbidden [O III] emission from the WIM, which show x(0 +2 ) < 0.1 in all cases and < 0.05 for 
a majority of the diffuse gas (Madsen et al. 2006). O III probes a range of ionization energies 
35.12 to 54.93 eV, similar to those probed by S +3 (see Figure[5]). Similarly, observations of He I 
recombination radiation limit the ionization fraction of singly-ionized helium. Observations find 
ionization fractions of singly-ionized helium as high as x(He + ) < 0.6, although many sight lines 
have x(He + ) < 0.3 or so (Madsen et al. 2006; Reynolds & Tufte 1995). Given the lack of higher 
ionization states of S from our results, the remainder of the He is likely to be neutral (see £J5]). 

There seems to be a significant amount of variation in the +2 and He + ionization fractions 
with location in the Galaxy. The highest ionization fractions from Madsen et al. (2006) in particular 



5 While we do not have a probe of S +4 , the lack of significant quantities of any other species with ionization 
energies > 35 eV (in the WIM or the hot gas up to > 150 eV; see below), we do not expect S +4 to be present in 
significant amounts. 
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seem to be found associated with energetic regions very near the disk or in unusual filaments. The 
diffuse WEVl and even supershells tend to trace the lower end of these ionization fractions. Figure 
[5] shows the full range of observed values, from the lowest upper limits (shown with downward 
pointing arrows) to the more extreme, higher ionization sight lines (Madsen et al. 2006). However, 
on the whole the emission line observations paint a similar picture to our absorption line study: 
ions requiring > 35 eV for production make a relatively small contribution to the total. 

While we have concentrated to this point on the ionization fractions within the WIM, we also 
have an opportunity along these sight lines to consider the ionization fractions of a broad range of 
ions in the total ionized gas of the Milky Way (i.e., WIM+HIM). Here we write the total column 
density of H II as 

iV(H II) = [(1 - f H IM)VWIM + fHIMVHIAl]N(e~), (6) 

where /him = iV(H II) him/N(H II) is the fraction of the H II column associated with the HIM. 
This value is derived from our estimates of the H II column associated with the hot ISM as dis- 
cussed in §2.11 While the HIM fraction is model-dependent and poorly constrained observationally, 
it is only important in this case so far as it affects the mean value of 77, which is given in the term in 
square brackets in Equation^ We assume 50% uncertainties in iV(H U)him, although its precise 
value does not affect the total H II column substantially. With this total H II column in hand, we 
write the total ionization fraction of an ion X J as 

x(x% tal = ^^a(X)-\ (V) 

equivalent to Equation 2. Figure [6] shows the results of this approach as applied to the S ions as 
well as the high ions Si +3 , C +3 , N +4 , +5 . The values are summarized in Table|4] We adopt solar 
system abundances for these ions. For CNO, these abundances are probably reasonable, since those 
elements are not depleted by more than a factor of 2 (perhaps less in the hot ISM). Incorporation 
of Si into the dust phase could lower the gas-phase abundance of Si quite a bit, probably even in 
the WIM (Howk & Savage 1999). Thus, the ion fraction of Si +3 could be somewhat higher than 
shown, perhaps by up to a factor of a few. 



5. Discussion 

We have presented direct measures of the ionization fractions of several ions of S in the 
WIM of the Milky Way (Figure HI), as well as of the ionization fractions of several metal ions in 
the integrated (warm and hot) ionized gas of the Milky Way (Figure [6]). These figures represent 
"maps" of the preferred ionization energies in the ionized gas of the Milky Way, showing at what 
energies we expect to find significant ionization fractions, albeit maps that are slightly skewed by 
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atomic physics (through varying ionization cross sections and recombination coefficients). In both 
cases, we find only those ions requiring < 35 eV for their production are present at more than the 
few percent level up to energies ~ 150 eV. This constrains the sources of ionization for the WEV1 
and the relative mass contributions as a function of temperature for the hotter, collisionally ionized 
gas (see below). 

The WIM shows no absorption from S ions above S +2 . The limits on S +3 and higher absorp- 
tion are meaningful for several reasons. First, the lack of high ionization S absorption suggests 
that the only ionization states of S with significant ionization fractions in the WIM are S + and 
S +2 , allowing us to estimate the S + ionization fraction (see $4] and Table [3]). The S + ionization 
fractions in the WIM along these two sight lines is quite a bit higher than found in classical O-star 
H II regions, which typically have a:(S + ) ~ 0.25 (Haffner et al. 2009), with S +2 being a dominant 
ion state in those regions. Thus, the WEVI has a generally lower ionization state than these H II 
regions, as noted previously (Haffner et al. 1999, Madsen et al. 2006). 

Second, the lack of S +3 (and S +5 toward M 3-vZ 1128) severely limits the amount of high 
ionization gas produced by the ionization sources for the WIM. The energies required to produce 
and destroy S + are 10.36 eV and 23.33 eV; for S+ 2 these are 23.33 eV and 34.83 eV. Thus, there is 
not a significant contribution to the S budget from any ions requiring > 35 eV for their production. 
Essentially all of the S in the WIM is in either S + or S +2 . While the ultimate sources of WEVI 
ionization are not precisely known, our measurements require they not produce significant amounts 
of S ions higher than S +2 . This is of particular importance for understanding the ionization state of 
He. The energy required to ionize S + is close to that required to ionize He (23.33 eV versus 24.59 
eV), and the ionization potential of S +2 is well below that required to ionize He + (47.22 eV versus 
54.42 eV). Thus, the lack of S +3 (probing ionization energies 34.79 to 47.22 eV) thus argues that 
virutally no He +2 (probing energies 54.42 eV and above) can be present in the WIM. 

This conclusion about the lack of He +2 is at odds with the work of Arabadjis & Bregman 
(1999), who claim that "there is little room for warm ionized gas of moderate ionization state." 
Based on an analysis of the ISM opacity to X-ray emission from distant galaxy clusters, these 
authors conclude that the fraction of neutral He in the WIM must be very small, with all of the He 
in either He + or He +2 . Since the emission line observations of He I recombination emission only 
provided a measure of the He + ion fraction, this conclusion could not immediately be ruled out by 
the existing observations. However, our results for the ionization fractions of S +3 and S +5 coupled 
with measurements of [O III] emission (Madsen et al. 2006, Reynolds 1985b) show that ions 
requiring > 35 eV for creation are trace ions in the WIM. Our limits on S IV give x(S +3 ) < 0.04 
(3<t) for both sight lines. This ion is destroyed at an energy of ~ 47 eV, well below the energy 
required to produce He + , ~ 55 eV. The ionization cross sections and recombination coefficients of 
S and He are not so different that a very high ionization fraction of He +2 could coexist with such 
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low ionization states of S +3 (e.g., Sembach et al. 2000). 

Several models have been developed in an attempt to match the metal ion emission from the 
WIM. Most of these are ID models relying on photons from OB stars to provide the ionization 
(e.g., Mathis 1986, 2000; Domgorgen & Mathis 1994, Sembach et al. 2000, Elwert & Dettmar 
2005). These are necessarily simplified models that often assume a single temperature stellar at- 
mosphere and very basic geometry, but include a great deal of atomic physics. Because many of 
these are largely concerned with matching the emission lines, which are temperature dependent, 
they do not all discuss directly the metal ion fractions in a way that is appropriate for comparison 
with the absorption line-derived results presented in this work. Sembach et al. (2000) attempted to 
model the WIM ionization using a series of ID models within Cloudy. They specifically tailored 
their models to assess the impact of the WIM on absorption line measurements. Their recom- 
mended composite model gives x(S + ) = 0.81 and x(S +2 ) = 0.18. The S +2 fraction, in particular, 
is lower by factors of 2 to 3 than our observations imply. Indeed, Howk et al. (2006) noted that 
applying the Sembach et al. (2000) models to the M 3-vZ 1 128 sight line did not produce results 
consistent with the full abundance determination for sulfur (underestimating the total S abundance 
by -0.4 dex). The sign of this result implies the ratio x(S + ) /x(S +2 ) in the Sembach et al. models 
is too large, consistent with our determinations. 

More recent models have considered the modification of the radiation field as it propagates 
through the ISM (e.g., Wood & Mathis 2004, Giammanco et al. 2004). The absorption of some 
radiation as it passes through the ISM will remove photons just above 1 Rydberg, "hardening" 
the H-ionizing photons as E ps 1 Ryd photons are removed while higher energy photons see less 
opacity. However, this will tend to increase the S + ion fractions compared with S +2 , since the 
former is created by unabsorbed E < 1 Ryd photons while the radiation field is diminished at 
the energy required to create S +2 (Hoopes & Walterbos 2003). However, this is true only when 
keeping all else constant. Elwert & Dettmar (2005) attempted to model the run of [S II] emission 
with height above the Perseus Arm observed by Haffner et al. (1999), following the propagation 
of radiation through the ISM with height. However, while their models do produce significantly 
higher S + ion fractions well above the plane (and hence are able to match the run of [S II] emission 
at heights z > 1 kpc above the Perseus Arm), they produce much lower S + ion fractions at low 
heights. Their models have x(S + ) = 0.2 to 0.5, with the ion fraction increasing with height up 
to z ~ 2 kpc (see their Figure 5). Integrating their S + ionization fractions vertically through an 
assumed exponential WIM distribution with a 1 kpc scale height (e.g., Haffner et al. 1999) gives a 
mean (x(S + )} = 0.28, lower by a factor of 2 than the value implied by our observations. A lower 
latitude sight line through such a model (approximating that to M 5-ZNG 1) would yield an even 
lower mean ionization fraction. The majority of the S in these models is likely in S +2 at z < 1 kpc, 
which dominates the column density of the distribution. The S +2 ion fractions in these models are 
significantly higher than our measurements imply. 
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We are not aware of significant constraints on the S +2 ionization fraction from emission line 
observations that could have been used to guide these earlier models. The emission line diagnostics 
also have an additional dependence on temperature that can mask mismatches between the model 
and WIM ionization fractions. Our measurements provide complementary constraints for models 
of the WIM. The causes of the discrepancies between the two models discussed above and the 
ionization fractions measured here are unknown at this point, but may include inappropriate choice 
of ionizing spectra, effects associated with the 3D structure of the gas and radiative transfer, and 
the lack of appropriate heating sources (which affects the match with emission line diagnostics, 
but not the absorption line diagnostics here). The ionizing spectrum is likely to be quite complex 
with multiple types of sources contributing. While it is almost certainly dominated by O star 
radiation (Haffner et al. 2009), that radiation is processed by intervening absorption and there may 
be contributions from other sources such as stellar remnants (e.g., Bregman & Harrington 1986, 
Rand et al. 201 1) and cooling radiation from hot gas and transition-temperature interfaces (Slavin 
et al. 2000). 

We note that the conditions in the WIM near the sun may be somewhat different than those 
seen in other galaxies, where obserations of the diffuse ionized gas (DIG) often trace layers brighter 
than the WIM of the Milky Way (Haffner et al. 2009). In particular, Rand et al. (2011) have 
recently used Spitzer observations of the [Ne III]/[Ne II] ratio to show the ratio of Ne +2 /Ne + is 
increasing with height above the plane. Because these emission lines are not very sensitive to 
excitation in the way the optical forbidden lines can be, this is a robust measure of the relative 
importance of these Ne ions. Observations of [O III] emission also find some galaxies show higher 
+2 ionization fractions than implied by the WHAM measurements of the Milky Way. Collins & 
Rand (2001) study, for example, showed that +2 could be the dominant ionization state of O in 
the extraplanar DIG of NGC 5775 and UGC 10288. One concern with comparing the Milky Way 
and these other galaxies is the difference in the properties of the observed DIG layers compared 
with the Milky Way's WIM. The former tend to be quite bright compared with the WIM and are 
observed at projected radii well within the solar circle in most cases. However, this is suggestive 
that the results found here may not apply across all galaxies and even over all positions within the 
Milky Way. 

The ionization fractions shown in Figure [6] demonstrate the importance of the WIM in the 
strong peak between ~ 10 and 35 eV, but they also shed light on the conditions in the hotter ionized 
gas of the Milky Way. With the exception of S + and S +2 , none of the ions in Figure [6] exceeds an 
ionization fraction of 0.05. The very low ionization fractions of the high ions Si +3 , C +3 , N +4 , +5 
have two causes. First, none of these Li-like ions (or Na-like in the case of Si IV) is expected to be 
the dominant ionization state at any temperature, having peak ionization fractions of rj 0.25 (Gnat 
& Sternberg 2007). However, the values seen in Figure [6] are an order of magnitude or more lower 
than this. The low ionization fractions for these high ions reflect the relatively small contribution 
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of gas in the transition temperature regime with T ~ 7 x 10 4 — 4 x 10 5 K to the total ionized gas 
content of the Milky Way. This is understandable: gas at these temperatures is at the peak of the 
interstellar cooling curve and radiates its energy very quickly, thereby cooling very quickly through 
this temperature regime. Thus, while the total mass budget of such transition temperature gas is 
not large compared with the WIM or the hotter HIM (see below), it is very important energetically. 
The mass flux or cooling flux through this regime is likely to be quite important. In fact, the ion 
fractions for these high ions include contributions from such transition-temperature gas as well as 
the ~ 40% of such gas that is at T < 7 x 10 4 K, representing matter that has already cooled and is 
out of ionization equilibrium or gas that is photoionized by the radiatively cooling HIM (Lehner et 
al. 2011). 

Figure [6] hints that there are likely to be two peaks in the total ionization fractions in the 
ionized gas of the Milky Way. One is associated with the WIM and traced by ions that probe the 
energy range E ~ 10 to 35 eV, as demonstrated by the large ionization fractions of S + and S +2 . 
This gas is predominantly photoionized and warm. The second peak is unseen in this plot, since we 
have no ions probing energies E > 140 eV. The gas associated with the hot ISM at T > 4 x 10 5 K 
(above that typically probed by O VI) will be traced in O by +6 and +7 ; such gas will be traced 
in S by S +6 and higher (e.g., Sutherland & Dopita 1993, Gnat & Sternberg 2007). While we cannot 
directly probe these ionization states along these sight lines, we can strongly limit their combined 
contribution. Summing the ionization fractions of S + and S +2 , we have x(S + ) + a;(S +2 ) ~ 0.80 ± 
0.10 along both the M 3 and M 5 sight lines. Thus, the combined fractions of all other ionization 
states represent only « (20 ± 10)% of the total S in these directions. The second peak in the 
ionization fractions of S ions should trace the high temperature HIM at a few hundred eV and have 
a magnitude of < 0.2. Figure [7] shows a plot of our derived ionization fractions supplemented by 
a calculation of sulfur ionization fractions for a hot, collisionally-ionized phase with T = 2 x 10 6 
K (e.g., Yao & Wang 2006). The hot phase ion fractions are from a Cloudy (Ferland et al. 1998) 
calculation assuming collisional ionization equilibrium (CIE) and have been normalized to sum to 
20% of the total (although the total for these high energies may be somewhat less). Thus, this figure 
shows a map of the expected energies at which one expects to find ions within the ionized medium 
of the Milky Way. The distribution in the highest ions depends on the assumed temperature of the 
HIM, but the other ions now have firmly measured ionization fractions along these two sight lines. 

The very low ion fractions of the high ions observable in the UV (esp. +5 ) have implications 
beyond the Galaxy. The strong O VI doublet is among the best tracers of highly-ionized gas 
in galaxy halos over all redshifts (Prochaska et al. 2011, Tumlinson et al. 2011). Studies of 
this sort that attempt to make estimates of the mass of gaseous galaxy halos seen in absorption, 
correcting the O VI columns based upon an assumed ionization fraction from collisional ionization 
models (e.g., Gnat & Sternberg 2007). In these cases, the typical ionization fraction adopted is 
a;(0 +5 ) < 0.2 (Gnat & Sternberg 2007, Sutherland & Dopita 1993). This may be the ionization 
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fraction in the O Vl-bearing gas at the temperatures where O VI peaks in abundance, and this gas 
represents a small fraction of the total, as we have discussed above. However, if one wishes to 
use O VI to calculate the total mass of ionized gas or the mass of hot ionized gas (e.g., T > 10 5 
K in this case), our results suggest some caution may need to be used, asa much lower mean 
ionization fraction is suggested for our two sight lines. In the two sight lines studied here, we 
find x(0 +b ) < 0.01 when considering the warm and hot ionized gas together. The O VI along 
these sight lines represents < 5% of the oxygen in the hot gas. The applicability of our results, 
which probe the thick disk of the Milky Way within z < 10 kpc of the midplane, to much different 
environments is unclear. Many studies in which such corrections are important probe galaxy halos 
at several tens of kpc impact parameter from galaxies (e.g., Prochaska et al. 2011, Tumlinson et 
al. 201 1) or probe the outskirts of HVCs (e.g., Fox et al. 2006, Sembach et al. 2003), which may 
have different ionization conditions than the gas in the Milky Way thick disk. 



6. Summary 

We have presented a method for studying the ion fractions of metal ions in the ionized gas of 
the Galaxy, expanding upon the method presented by Howk et al. (2006) for measuring elemental 
abundances. We make use of unique sight lines toward globular clusters containing both UV bright 
stars (for measuring the metal ions) and pulsars (for measuring the electron column density). We 
apply this method to study the ionization of the ionized gas toward M 3-vZ 1 128 and M 5-ZNG 1, 
both of which probe extended paths through the Galactic thick disk or halo, providing estimates 
for the ionization fractions of S°, S + , S +2 , S +3 , and S +5 in the WIM. We also assess the ionization 
fractions integrated through the warm and hot ISM of the Milky Way. 

Our principal conclusions based on this analysis are as follows. 

1. The only ions with significant ionization fractions in the WIM are those with ionization 
potentials < 35 eV. We find S +2 makes up a substantial portion of the WIM sulfur, with 
ionization fractions x(S +2 ) = 0.33 ± 0.07 and 0.47 ± 0.09 toward M 3 and M 5, respectively. 

2. We limit the contribution of S +3 and higher ionization states of sulfur, with 3a upper limits 
x(S +3 ) < 0.06 for both of our sight lines. Given the lack of higher ionization states in the 
WIM, the ionization fractions of S + and S +2 represent the vast majority of the WIM sulfur. 
Assuming these add to unity, we derive x(S + ) = 0.67 ± 0.07 and 0.53 ± 0.09 toward M 3 
and M 5, respectively. These numbers are in good agreement with emission-line derived S + 
ion fractions (Haffner et al. 1999, Madsen et al. 2006). 

3. Existing, simplified photoionization models for the WIM produce S ion fractions that can 
differ by factors of 2 from our estimates, with various models under- or over-producing S +2 . 
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These models are largely constructed to match emission line measurements, which provide 
access to a narrower range of ions than our measurements. The lack of higher ionization 
states in the WIM rules out the conclusion of Arabadjis & Bregman (1999) that much of 
the He in the WIM is in the form of He +2 and limits the hardness of the radiation field 
responsible for ionizing the WIM. 

4. We show that the ionized gas column density (and hence mass) along these two sight lines is 
dominated by warm photoionized gas (T ~ 10 4 K) favoring ionization states with ionization 
energies E ~ 10 to 35 eV and a hot ionized phase (T > 4 x 10 5 K) favoring energies of a 
few hundred eV. The warm phase contains ~ 80% of the total ionized gas column, and the 
hot phase < 20%. The important tracer O VI probes only a small fraction of the ionized gas 
along the sight line with x{0 +b ) < 0.01. 
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Fig. 1.— Absorption line profiles of S II AA 1250.584, 1253.811, 1259.519, and S III A 1190.208 
from STIS E140M observations as well as S III A 1012.495 from FUSE observations of M5-ZNG1 . 
The lighter colored portion of the spectra are interloping lines. The FUSE data are shown with a 
binned pixel size of 3.84 km s _1 , giving ~ 5 pixels per resolution element. The STIS data are 
shown with their native pixels, that is 3.22 km s -1 , or ~ 2.3 pixels per resolution element at these 
wavelengths (Proffitt et al. 2002). 
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Fig. 2. — Comparison of the apparent column density profiles for S II A1250.584, 1253.811, and 
S III A1190.208 from STIS E140M observations. The top panel is a comparison of the two S II 
profiles. The discrepancy between the two profiles is due to the presence of unresolved saturation. 
The lower panel compares the S II line at 1250.54 A with the S III line at 1190.208 A. The S III 
profile has been "cleaned" of Si II contamination. 
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Fig. 3. — Two views of the STIS E140M spectrum of the Lyman-a; absorption profile toward M 5- 
ZNG 1. Both stellar and interstellar absorption contribute to this profile. The estimated stellar 
profile is shown as the thin blue line and has been shifted to t> L sR = +65.7 km s _1 to match the 
observed positions of stellar absorption lines in the STIS spectrum. The best fit to the interstel- 
lar and stellar absorption profile is shown in red, corresponding to an interstellar column density 
logiV(H I) = 20.47 ± 0.02. The sharp line in the center of the Lyman-a absorption trough is 
geocoronal emission. 
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Fig. 4. — The ionization fractions, x(X^), of S ions in the WIM as a function of ionization energies 
for the M 3 (black) and M 5 (red) sight lines. The horizontal extent of the bars for each ion span 
the range of energies required to create it, from the next lower ionization state, and destroy it, 
creating the next higher ionization state (i.e., the ionization energies). The M 5 ionization energies 
are offset by 0.5 eV for clarity. 
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Fig. 5. — A comparison of emission and absorption line constraints on ionization fractions for 
several ions in the WIM. The ionization fractions derived here for S ions are shown for the M 3 
(black) and M 5 (red) sight lines as in Figure 0] The shaded regions denote the range of values 
seen by emission line constraints on the ionization fractions of of S + (Haffner et al. 1999, Madsen 
et al. 2006 as summarized in Haffner et al. 2009), +2 (Madsen et al. 2006, Reynolds & Tufte 
1995), and He+ (Madsen et al. 2006, Madsen 2004, Reynolds 1985b). We have excluded the 
values associated with the northern filament studied by, e.g., Madsen et al. (2006) that seems to 
show higher ionization fractions for +2 and He + than typically found in the diffuse WIM. The 
lowest upper limits on the emission constraints are shown with downward facing arrows. 
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Fig. 6. — The ionization fractions integrated over all of the ionized gas (WIM+HEVl) of the Milky 
Way for several metal ions as a function of ionization energies for the M 3 (black) and M 5 (red) 
sight lines. This plot differs from Figure|4]in that it applies to both the WIM and HIM. The results 
are tabulated in Table|4l The mass of ionized gas (at least along these sight lines) is dominated by 
WIM gas traced by ions in the range ~ 10 to 35 eV and by a hotter component with T > 4 x 10 5 
K probed by ions with higher ionization energies than probed by the available UV transitions. To- 
gether the WIM ions account for 80% of the total, implying that ions tracing the high temperature 
phase cannot represent more than ~ 20% of the total contribution to any metal. Gas with temper- 
atures in the range a few x 10 4 K to ~ 4 x 10 5 K contributes very little to the total mass budget, as 
discussed in §5\ 
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Fig. 7. — The ionization fractions for the total ionized gas (WIM+HIM) of the Milky Way, as 
Figure [6l for several metal ions as a function of ionization energies (plotted on a log scale here), 
including CIE models for S ions in the hot ionized phase (blue). The CIE models assume T = 
2 x 10 6 K (e.g., Yao & Wang 2006) for the HIM and have been normalized to contribute (20 ±5)% 
of the total ionization fraction, although the contribution could be somewhat lower. The ionization 
fractions of all the S ions in this figure sums to unity. This plot represents a map of the energies at 
which we expect significant ionization fraction contributions over all species. The distribution of 
S ions in CIE will, of course, depend on the assumed temperature. 
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Table 1 . Adopted Column Densities Toward M 3 and M 5 





M3-vZ1128 


M 5 - ZNG 1 


Species 


logiV Ref. a 


logiV Ref. a 



HI 


19.98 ±0.03 


1 


20.57 ±0.03 


4 


e~ 


19.91 ±0.01 


2 


19.96 ±0.01 


2,5 


C IV 


14.39 ±0.03 


4 


14.38 ±0.02 


4 


N V 


i q op;+0.07 
io - /d -0.09 


4 


13 in+ - 10 

lO.±U_ .l 4 


4 


VI 


14.49 ± 0.03 


3 


14.41 ±0.02 


4 


Si IV 


13.80 ±0.02 


4 


13.79 ±0.02 


4 


S I 


< 12.7 (3a) 


3 


< 12.7 (3a) 


4 


S II 


15.28 ±0.02 


3 


15.59 ±0.10 


4 


S III 


14.47 ±0.03 


3 


14.66 ±0.02 


4 


S IV 


< 13.7 (3a) 


3 


< 13.7 (3a) 


4 


S VI 


< 13.4 (3a) 


3 







References: (1) Howk et al. 2006 ; (2) Hessels et al. 
2007; (3) Howk et al. 2003; (4) This Work; (5) Freire et al. 
2008. 
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Table 2. Derived Interstellar Parameters Toward M 3 and M 5 



Quantity M 3 M 5 



logiV(e - ) 


19.91 ±0.01 


19.96 ±0.01 


log iV(H \\) HIM 


19.20 ±0.18 


18.93 ±0.18 


log N(R R) WIM 


19.79 ±0.06 


19.83 ±0.06 


logiV(H) a 


20.24 ±0.03 


20.66 ±0.02 


iV(H H)/JV(H) 
logA(S) b 


0.45 ±0.09 


0.19 ±0.03 


-4.86 ± 0.04 


-5.00 ±0.10 


[S/H] 


-0.01 ±0.04 


-0.15 ±0.10 



a The total hydrogen column density, including 
contributions from both warm and hot gas, i.e., in- 
cluding the ionized gas associated with the HIM. 

b The sulfur abundance is derived excluding the 
H II column from the HIM. Thus, the hydrogen ref- 
erence columns for comparison with the summed S 
ion column densities are logiV(H) = 20.20 ± 0.03 
and 20.64 ± 0.02 for M 3 and M 5, respectively. 
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Table 3. WIM Sulfur Ionization Fractions a 



Quantity 


M3 


M5 b 


x(S°) 


< 0.006 (3a) 


< 0.005 (3a) 


x(S + ) 


0.67±0.07 c 


0.53±0.09 c 


x(S +2 ) 


0.33 ±0.07 


0.47 ±0.09 


x(S +3 ) 


< 0.06 (3a) 


< 0.05 (3a) 


x(S+ 5 ) 


< 0.03 (3a) 





a The ionization fractions are x(X^) = 
N(Xi)/N(X). We assume x(R + ) = 
0.95 ± 0.05 in deriving these values (see 
text). 

b We adopt the abundance towards M 3 
for the ionization fraction calculations. If 
one adopts the less well-determined value 
derived for M 5, the result for S +2 , for ex- 
ample, is x(S +2 ) = 0.67 ± 0.19. 

c The S + ion fractions are not measured 
directly; instead they are derived assuming 

x(S + ) + x(S +2 ) = 1 (see text). 



Table 4. Total (WIM+HIM) Ionization Fractions 



Quantity M 3 M 5 



s(S°) 


< 0.005 (3<r)° 


< 0.004 (3a) c 
0.42±0.10 d 


x(S+) 


0.53±0.10 d 


x(S+ 2 ) 


0.27±0.04 c 


0.38±0.08 c 


x(S+ 3 ) 


< 0.05 (3(t) c 


< 0.04 (3(j) c 


x(S+ 5 ) 


< 0.02 (3(j) d 




x(C+ 3 ) 


n m O+0.002 d 

U.U±Z,_Q QQQg 


n in+0.001 d 
u - lu -0.0007 


x(N+ 4 ) 


0.003 ±0.001 e 


0.002 ± 0.002 e 


x(0+ 5 ) 


0.008 ±0.001° 


0.006 ±0.001° 


x(Si+ 3 ) 


0.025 ±0.03° 


0.022 ± 0.003° 



a The ionization fractions are x(X j ) = 
N(X')/N(X). We assume x(R+) = 0.95 ± 
0.05 (see text). For the hottest gas, we expect 
this fraction to be unity. 

b We adopt the abundance towards M 3 for 
the ionization fraction calculations. 

d The S + ion fractions are not measured 
directly; instead they are derived assuming 

x(S + ) + x(S +2 ) = 1 (see text). 

°Ions higher than S +3 have an assumed 

x(H+) = 1.0 



